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resultado de las practicas forestales y agricolas, y se espera que ambas 
cinco impactos diversos en los servicios del ecosistema (ES). En este 
estudio, investigamos esta declaracion general para regular y apoyar los 
servicios a traves de una evaluacion de las funciones del ecosistema: 
regulacion climatica a traves del secuestro de carbono en el suelo y la 
biomasa vegetal, ciclo del agua y regulacion de la erosion del suelo a 
traves de la infiltracion de agua en el suelo, y apoyo a la produccion 
primaria a traves de calidad qulmica del suelo y almacenamiento de agua. 
Probamos la hipotesis de que los patrones de composicion y estructura de 
la cobertura del suelo alteran significativamente las metricas de ES en 
dos escalas diferentes. Examinamos 54 granjas en dos regiones amazonicas 
de Brasil y Colombia y evaluamos la composicion y estructura de la 
cubierta terrestre a partir de dates de deteccion remota (escala de 
granjas) de 1990 a 2007. Se utilizaron metodos simples y bien 
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(2014) 29: 311-328 DOI 10.1007 / sl0980-013-9981-y desde cinco puntos en 
cada granja (escala de parcela). La mayoria de las metricas de ES se 
correlacionaron significativamente con las clasificaciones de uso de la 
tierra (escala de parcela) y de cobertura de la tierra (escala de 
granja); sin embargo, la variabilidad espacial en las propiedades 
inherentes del suelo, sola o en interaccion con el uso del suelo o los 
Cambios en la cubierta del suelo, contribuyo en gran medida a la 
variabilidad en las metricas de ES. Las existencias de carbono en la 
biomasa de plantas por encima del suelo y la tasa de infiltracion de agua 
disminuyeron de bosques a pastizales, mientras que la calidad quimica del 
suelo y la capacidad de almacenamiento de agua disponible para las 



plantas aumentaron. Las clasificaciones de cobertura del suelo basadas en 
metricas de estructura explicaron significativamente menos variacion de 
metrica ES que aquellas basadas en metricas de composicion. La dinamica 
de la composicion de la cobertura del suelo explico el 45% (P \ 0.001) de 
la varianza metrica ES, el 15% por si mismo y el 30% en interaccion con 
las propiedades inherentes del suelo. Este estudio describe como 
evolucionan los ES con los cambios del paisaje, especificando la 
contribucion de la variabilidad espacial en el entorno fisico y 
destacando las compensaciones y las sinergias entre los ES. Palabras 
clave Intensidad del uso de la tierra Servicios del ecosistema del suelo 
Impulsores socioeconomicos Agroecosistemas Almacenamiento de carbono 
Calidad quimica del suelo Infiltracion del agua Compensacion Introduccion 
La Evaluacion de Ecosistemas del Milenio (2005) llamo la atencion sobre 
la degradacion y el uso no sostenible de los servicios del ecosistema 
(ES) esenciales para la vida humana. Comprender los procesos ecologicos 
detras de ES y medirlos aun plantea muchas preguntas cientificas 
(Carpenter et al. 2006; Bennett et al. 2009) . Nuestra capacidad para 
disenar paisajes que suministren multiples ES de manera sostenible 
adolece de una falta de conocimiento sobre los procesos ecologicos detras 
de los ES y sus respuestas al uso de la tierra y los cambios globales 
(Foley et al. 2005; Carpenter et al. 2009; Turner 2010) . Los impulsores 
socioeconomicos y ecologicos involucrados en el suministro de ES, las 
relaciones entre ES y los mecanismos subyacentes no se conocen bien 
(Kremen y Ostfeld 2005; Nelson et al. 2009; Breure et al. 2012) . Una 
mejor comprension parece critica para evitar intercambios no deseados 
entre diferentes ES y, a la inversa, para beneficiarse de sinergias 
positivas entre ellos (Bennett et al. 2009). La region amazonica es un 
lugar ejemplar para abordar estas preguntas (Foley et al. 2007). Muy 
diversa (Hoorn et al. 2010), el bosque amazonico tambien se considera una 
de las ultimas reservas de tierra para satisfacer las crecientes 
necesidades humanas de alimentos y otros productos agroforestales (Foley 
et al. 2007) . Aunque la agricultura no es el unico impulsor directs del 
desarrollo regional, sigue siendo la principal causa de deforestacion 
(Fearnside 2005; Aldrich et al. 2006; Godar et al. 2012) . La agricultura 
de tala y quema practicada en el frente pionero se reemplaza 
progresivamente por pastos y cultivos de exportacion en etapas 
posteriores (Rodrigues et al. 2009) . Sin embargo, los patrones y la 
dinamica del uso de la tierra varian en gran medida, dependiendo 
principalmente de las condiciones socioeconomicas que prevalecen a nivel 
local y regional. Incluso la deforestacion parcial por la tala o el fuego 
representa una amenaza inmediata para la biodiversidad (Gibson et al. 
2011) y el clima (Malhi et al. 2008) . El ecosistema forestal es 
resistente (es decir, puede mantener su composicion y funciones) cuando 
el deficit de agua estacional es moderado (Phillips et al. 2009). Sin 
embargo, debido a las interacciones entre la deforestacion, los incendios 
forestales y la sequia, no se puede descartar que el sur y el este de la 
cuenca del rio Amazonas puedan estar en un estado transitorio hacia un 
sistema perturbado, propenso a la liberacion de carbono, que afectaria el 
agua y los productos biogeoquimicos. ciclos (Davidson et al. 2012). 
Estrechamente relacionada con la perdida de biodiversidad, la degradacion 
del suelo en los pastos no es menos preocupante (Mathieu et al. 2005, 
2009), con consecuencias para los ciclos de agua y nutrientes y la 
produccion primaria (Markewitz et al. 2004; Martinez y Zinck 2004; Chaves 
et al.2008). Enfrentados a una disminucion en la productividad de los 
pastizales, los pequenos propietarios generalmente abandonan las tierras 
que limpiaron para granjas grandes y finalmente continuan su migracion, 
contribuyendo asi al avance del frente de deforestacion D. Unidad de 
Investigacion del Espacio para el Desarrollo D. Mitja (ESPACE-DEV UMR), 
Instituto de Research for Development (IRD), 34093 Montpellier, Francia 
M. Sarrazin Institute of Research for Development (IRD), 97323 Cayenne, 
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actores buscan innovaciones tecnicas o reclaman nuevos incentivos de 
polltica publica, como los mecanismos de compensacion (por ejemplo, '' 
creditos de carbono '') recomendados en las negociaciones para reducir 
las emisiones de la deforestacion y la degradacion forestal (REDD) (Hall 
2008). El future de la region amazonica depende de la capacidad de los 
ecosistemas para soportar las perturbaciones causadas por el use de la 
tierra y los cambios climaticos y en la implementacion de pollticas para 
conciliar el desarrollo y la conservacion (Betts et al. 2008) . Para 
promover la resiliencia de los ecosistemas influenciados por el ser 
humane y mantener los beneficios globales proporcionados por el bioma 
amazonico, se deben brindar servicios de apoyo y regulacion, 
proporcionados por las funciones del ecosistema, como la produccion 
primaria, el reciclaje de nutrientes, el secuestro de carbono y la 
infiltracion de agua en el suelo. sostenido. En este artlculo, analizamos 
como estos servicios, evaluados a traves de funciones, evolucionan con 
los cambios del paisaje a lo largo de un gradiente de intensidad de use 
de la tierra. Probamos la hipotesis de que los patrones de composicion y 
configuracion de la cobertura del suelo alteran significativamente la ES. 
Los ES abordados son la regulacion del clima, a traves del secuestro de C 
en el suelo y la biomasa vegetal, la regulacion del ciclo del agua y el 
control de la erosion, a traves de la infiltracion de agua en el suelo y 
el apoyo a la produccion primaria, a traves de la calidad qulmica del 
suelo y el almacenamiento de agua. Nuestros objetivos especlficos fueron 
(1) evaluar las contribuciones relativas de los sucesivos cambios en el 
uso de la tierra y las variaciones espaciales en el entorno flsico en los 
cambios de ES en la parcela y las escalas de la granja, (2) evaluar el ES 
en la escala de la granja utilizando metricas de tierra abarcan la 
composicion y / o configuracion y su dinamica, y (3) para revelar 
compensaciones y sinergias entre los ES durante la transformacion del 
paisaje y los procesos ecologicos en juego. Materiales y metodos Sitios 
de estudio El estudio se llevo a cabo en dos regiones amazonicas de 
Colombia (Departamento de Caqueta) y Brasil (Estado de Para) con 
historias de colonizacion distintas. Se seleccionaron tres areas 
representativas de la diversidad de los sistemas agricolas de pequenos 
productores en ambas regiones a lo largo de un gradiente de extension de 
la deforestacion y la intensidad del uso de la tierra agricola (Fig. 1). 
Queda poca cobertura forestal en las areas de estudio colombianas (20-60 
km de Florencia, region de piedemonte de la cordillera oriental andina) 
donde comenzo la deforestacion en la decada de 1940 (Tabla 1). Las tres 
areas diferlan en su manejo: los sistemas tradicionales en Canelos 
dedicados principalmente a la produccion de lacteos y la ganaderla 
extensiva en granjas relativamente grandes (Cuadro 1) con pasturas 
degradadas establecidas desde hace mucho tiempo; sistemas silvopastoriles 
en Balcanes basados ??en cultivos de leguminosas forrajeras para 
restaurar pasturas y mejorar la produccion de ganado; y sistemas 
agroforestales en Aguadulce que combinaban cacao, arboles frutales y 
pastos (Arnauld de Sartre et al. 2011). Las areas brasilenas (centre del 
estado de Para') fueron colonizadas mas recientemente y el 40-70% del 
bosque inicial permanecio en 2007 (Silva Costa et al. 2012) . La mayorla 
de los colonos del area de Mac.aranduba abandonaron sus actividades agro- 
extractivistas iniciales para dedicarse a la produccion ganadera debido 
al predominio de esta actividad y al agotamiento de los recursos 
madereros cerca de la ciudad de Maraba ', ambas fuentes de conflicto con 
los madereros y ganaderos vecinos. . En Palmares II, los agricultores 

obtienen ingresos de los cultivos alimentarios, en lugar de la produccion 
ganadera, para alimentar la cercana ciudad de Parauapebas. El area de 
Pacaja, un primer frente pionero que comienza desde un sendero 



(Travessa'o 338 Sul) perpendicular a la carretera TransAmazoniana, se 
dividio entre cultivos anuales y perennes (cacao) y produccion ganadera. 
Ademas, las dos regiones difieren notablemente en sus entornos flsicos. 
Las areas de Caqueta se encuentran en sedimentos del Mioceno con un 
paisaje relativamente piano, mientras que las areas de Para estan 
situadas en rocas cristalinas granlticas o metamorfleas de edad 
paleoproterozoica, con paisajes montanosos de aspereza media en 
elevaciones bajas (maximo 350 m en Brasil vs. 1,100 m en Colombia). Ambas 
regiones tienen un clima tropical humedo, pero la precipitacion anual 
promedio es de aproximadamente 1.700 mm en Para' y 3.300 mm en Caqueta. 
Diseno de muestreo Este estudio se basa en dates recopilados en 54 
granjas (nueve por area) elegidas de un conjunto inicial de 304 granjas 
(51 por area; faltan dates para dos granjas), a partir de las cuales se 
desarrollo una clasificacion de los sistemas agricolas, basada en una 
encuesta socioeconomica (Arnauld de Sartre et al. 2011) . Los grupos de 
sistemas agricolas se separaron mediante analisis jerarquico de 
conglomerados: granjas de Landscape Ecol (2014) 29: 311-328 313 123 Fig. 

1 Ubicacion de las areas de estudio y diseno de muestreo 314 Landscape 
Ecol (2014) 29: 311-328 123 pionero temprano los frentes se deforestaron 
hace 10-20 anos con una produccion mixta de cultivos y ganado anuales y 
perennes, y bajos ingresos; sistemas similares con ingresos ligeramente 
mejores; granjas con gran produccion ganadera y los mayores ingresos; y 
sistemas agroforestales establecidos en sitios donde la deforestacion 
comenzo hace 50-80 anos, con ingresos intermedios. Cada area se dividio 
en tres subareas, en las que se seleccionaron tres granjas no adyacentes 
para representar mejor la diversidad de los productos agricolas, la 
dinamica de la cobertura del suelo y el tamano de la granja identificados 
por encuestas socioeconomicas y de teledeteccion. Para cada una de las 54 
granjas, cinco puntos de muestreo estaban espaciados equitativamente a lo 
largo de un transecto correspondiente a la diagonal mas larga de la 
granja o a un eje mas o menos norte-sur (Fig. 1). La distancia entre 
puntos era, por lo tanto, igual a 1/6 de la longitud del transecto y 
variaba segun el area de la granja. Este enfoque aseguro que el esfuerzo 
de muestreo de los diferentes usos de la tierra fuera proporcional a su 
ocurrencia relativa en cada granja. Los puntos de muestreo fueron 
parcelas demarcadas de 50 9 10 m2 cuya longitud era perpendicular al 
transecto de muestreo de la granja. Descripcion de los usos del suelo a 
escala de parcela Registramos los usos del suelo en el campo al mismo 
tiempo que recolectamos dates para la evaluacion de ES. La combinacion 
posterior de algunas de las subcategorias identificadas en el campo (por 
ejemplo, tipos de cultivos perennes) permitio que cada una de las 11 
categorias de uso de la tierra utilizadas para el analisis estadistico 
contuviera al menos 10 observaciones (Tabla 2). Analisis de paisajes a 
escala de granja La mayoria de los metodos analiticos detallados estan 
disponibles en Oszwald et al. (2011) Conservamos a continuacion los 

aspectos necesarios para comprender las clasificaciones de cobertura del 
suelo utilizadas en nuestro estudio. Utilizamos datos de teledeteccion 
para caracterizar la dinamica de la cubierta terrestre de las 54 granjas. 
£,Usamos Landsat TM y ETM? imagenes (resolucion espacial de 30 m, 
grabacion espectral adaptada a la identificacion de la cobertura del 
suelo) tomadas durante las estaciones secas de 1990, 1994, 1998, 2002 y 
2007. Se aplicaron correcciones radiometricas y atmosfericas al conjunto 
de datos satelitales con el modelo 5S (Kergomard 2000) . Las mediciones de 
validacion de campo se tomaron durante la estacion seca de 2007 y 2008 
para construir un sistema supervisado de clasificacion de la cobertura 
del suelo con un algoritmo de maxima verosimilitud, que se vinculo a la 
firma espectral de cada una de las unidades de cobertura del suelo 
(Fig.2). La clasificacion supervisada se realize en la imagen Landsat TM 
de 2007. Las firmas espectrales permitieron la reconstitucion de las 
cubiertas de tierra pasadas de imagenes anteriores (1990, 1994, 1998 y 



2002). En cada fecha, la cobertura del suelo de cada granja se describio 
utilizando 10 metricas de composicion (proporciones relatives del area 
total de la granja Table 1 Ubicacion geografica y caracterlsticas 
principales de las areas de estudio Sitios X (UTM, WGS84) Min Max Y (UTM, 
WGS84) Min. Max. Inicio / llegada de la colonizacion 50% agricultores 
actuales (anos) Tamano promedio de la finca (ha) Cobertura forestal 
original restante en 2007 (%) Cobertura de pastures en 2007 (%) Colombia 
Canelos 2,311,340.8 10,168,671.4 1950/1991 184 3 66 2,320,397.3 
10,184,388.2 Aguadulce 2,324,742.8 10,187,231.1 1950/2000 49 2 39 
2,312,112.4 10,173,341.7 Balcanes 2,352,560.7 10,178,497.1 1940/2000 65 4 
72 2,335,720.1 10,163,170.7 Brasil Mac.aranduba 679,602.5 9,473,375.9 
1982/1992 81 39 17 686,282.9 9,467,677.1 Pacaja' (Trail 338S) 

493,189.79,519,2191919191919191919191919191919191919959199499,295,29,2199 
,219,292,292,2 II 625,263.4 9,343,195.9 1995/1997 25 38 18 630,403.4 
9,352,609.8 Paisaje Ecol (2014) 29: 3II-328 315 123 de las 10 unidades de 
cobertura terrestre presentadas en la Fig. 2) y nueve estructuras 
cumplidas rics, tanto de composicion como de configuracion (Lausch y 
Herzog 2002): area de cobertura de la tierra (area de la granja en ha), 
diversidad de parches (riqueza y diversidad de Shannon, indices de 
uniformidad y dominancia de parches), densidad media de parches (m ha-1) 

, densidad de borde (m ha-1) y estructuras fractales (relacion perimetro 
/ area, indice de forma media). Estas metricas se utilizaron para 
constituir matrices tridimensionales en las que las granjas se 
caracterizaron por su composicion y estructura de cobertura del suelo en 
las cinco fechas de las imagenes Landsat TM. El analisis de estas 
matrices con el analisis normalizado de componentes principales (PCA) y 
el metodo ACT-STATIS (Lavit et al. 1994) proporciono cuatro indicadores 
que se describen a continuacion (Fig. 2). El indicador de composicion 
claramente separo las fincas con predominio de bosque (valores cercanos a 
1) de las fincas con una alta proporcion de cobertura agricola (valores 
cercanos a 0). El indicador de estructura enfatizo las diferencias en la 
organizacion de la cobertura del suelo entre granjas, desde granjas en 
las cuales el mosaico de cubierta de tierra tenia una alta complejidad de 
forma en parches forestales (valores cercanos a 1) a otros que tenian 
parches agricolas grandes (valores cercanos a 0); los valores intermedios 
se asociaron con granjas que tienen alta complejidad en parches de 
vegetacion secundaria. El indicador de "dinamica de composicion" 
distinguia las granjas boscosas que habian experimentado poca o ninguna 
deforestacion desde 1990 (valores cercanos a 1) de las granjas con un 
paisaje agricola homogeneo desde al menos 1990 (valores cercanos a 0); 
los valores intermedios fueron una funcion de la intensidad (es decir, 
importancia y rapidez) de los cambios en la cubierta terrestre entre 1990 
y 2007. El indicador de "dinamica de la estructura" evaluo el cambio 
temporal en la organizacion de la cubierta terrestre y distinguio a las 
granjas que habian experimentado una deforestacion reciente con poca 
fragmentacion del bosque (valores cercanos a 1) de las granjas donde se 
produjo la deforestacion en la decada de 1990, seguida de una intensa 
transformacion del bosque a la agricultura (valores intermedios), y 
finalmente granjas con un paisaje agricola homogeneo (valores cercanos a 
0). Evaluacion del servicio del ecosistema Se estimaron cinco metricas 
(indicadores) de ES: (1) stock de C en biomasa de plantas (arbustos y 

arboles) y (2) stock de C del horizonte de suelo de 0-30 cm (servicio de 
regulacion climatica); (3) tasa de infiltracion de agua en el suelo 

(servicios de regulacion del ciclo del agua y la erosion del suelo); (4) 

un indice de calidad quimica del suelo y (5) capacidad de almacenamiento 
del suelo de agua disponible para la planta del horizonte del suelo de 0- 
10 cm (servicio de apoyo a la produccion primaria). Estas metricas de ES 
se estimaron en la escala de la parcela a partir de las caracteristicas 
de vegetacion y suelo medidas en todos los puntos (Tabla 3) y en la 
escala de la granja calculando los valores medios de los cinco puntos de 



muestreo dentro de una granja. Para las ultimas cuatro metricas, se 
excavaron cuatro pozos de suelo regularmente separados por 10 m en el 
medio de cada parcela de 50 9 10 m2. Las mediciones y la recopilacion de 
dates tuvieron lugar en 2008 durante 4 meses (abril-julio) de la 
temporada de Iluvias en cada pals. Se utilizaron metodos bien 
establecidos y relativamente simples, describes a continuacion, para 
caracterizar un total de 270 puntos. Existencias de C en la biomasa de 
plantas. Biomasa de arboles de plantas secas sobre el suelo [BT: diametro 
a la altura del pecho (dap) CIO cm] y arbustos (BB: dbh Tabla 2 Usos del 
suelo registrados en el campo en areas de estudio brasilenas y 
colombianas en el memento en que se obtuvieron los dates recolectada para 
la evaluacion metrica del servicio del ecosistema (ES) de un total de 270 
puntos (parcelas) y numero de puntos en cada clase de uso de la tierra y 
pals Clases de uso de la tierra Numero de puntos Brasil Colombia Basques 
preservados 15 (14) 0 Basques explotados 24 (21) 0 Basques quemados 10 0 
Cultivos anuales (arroz, yuca o malz) 14 0 Barbechos (basques secundarios 
en tierras de cultivo abandonadas) 17 0 Pastes limpiados 18 (17) 64 

Pastes invadidos 17 9 (8) Pastes mixtos (pastes con plantas nativas y 
sistemas silvopastoriles) 3 19 Cultivos perennes (cacao, plantaciones de 
arboles de rubro o sistemas agroforestales) 4 24 (22) Basques jovenes 
secundarios en pastures abandonadas 6 4 Basques viejos secundarios en 
pastures abandonadas 7 (6) 15 (14) Los numeros entre parentesis indican 

el numero de parcelas despues de eliminar 10 valores atlpicos puntos 
(valores extremes de metricas ES o propiedad inherente del suelo) 316 
Landscape Ecol (2014) 29: 311-328 123 \ 10 cm y altura [2 m) se estimo 
aplicando ecuaciones alometricas proporcionadas para los bosques (Higuchi 
et al. 1998; Gerwing 2002) y barbechos (Nelson et al. 1999) despues de 
identificar especies y medir el diametro y la altura de los individuos en 
parcelas de 50 9 10 m2 y 50 9 5 m2, respectivamente, y el contenido de 
agua a 70 C de las alicuotas de biomasa correspondientes ( todos los 
individuos con dap B 5 cm), de lo contrario Fig. 2 Unidades de cobertura 
terrestre distinguidas del analisis de teledeteccion (firma espectral 
Landsat TM, imagenes de 1990, 1994, 1998, 2002 y 2007) y validacion de 
campo (2007 y 2008) . Metrica del paisaje: cuatro indicadores del paisaje 
que describen la composicion y estructura de la cobertura terrestre en 
2007 y la dinamica de la composicion y estructura de la cubierta 
terrestre entre 1990 y 2007 Landscape Ecol (2014) 29: 311-328 317 123 
(dap [5 cm) aplicando el factor 0.603 segun Higuchi et Alabama. (1998) En 
ausencia de analisis de C de muestras de plantas, se uso un factor de 0,5 
para estimar la masa de C en la biomasa de plantas (Markewitz et al. 

2004). Stock de C en el suelo El contenido total de C del suelo en el 
horizonte de 0-30 cm se midio con un analizador CHNS en tres muestras de 
cada parcela (compuestos de las profundidades de 0-10, 10-20 y 20-30 cm 
de los cuatro pozos). Para calcular el stock de C, se repitieron las 
mediciones de densidades aparentes (qb; metodo de nucleo cilindrico, 7 cm 
de diametro) de 0 a 30 cm, separando profundidades de 0-2, 2-5, 5-10, 10- 
20 y 20-30 cm. en los cuatro hoyos Tasa de infiltracion de agua en el 
suelo Aplicando la prueba de infiltracion Beerkan (Labassate're et al. 
2006), calculamos la tasa de infiltracion (Vi) de un volumen fijo de agua 
de 250 cm3 vertido en un anillo simple de 20 cm de diametro insertado en 
el suelo superficie a una profundidad de aproximadamente 1 cm. Esta 
prueba de infiltracion se repitio cuatro veces, cerca de los pozos de 
tierra. Indice de calidad quimica del suelo La distribucion del tamano de 
particula y las propiedades quimicas del horizonte de suelo de 0-10 cm se 
determinaron en un compuesto de dos pozos. Las propiedades quimicas 
realizadas fueron pH H20, capacidad de intercambio cationico (CEC) a pH 
del suelo, intercambiable A13 ?, Ca2 ?, Mg2? y K ?, NH4 intercambiable? y 
fosforo extraible (metodo de extraccion de P, doble acido de Mehlich en 
HCl 0.05 M y soluciones de H2S04 0.0125 M) utilizando metodos estandar 
(Pansu y Gautheyrou 2006). Se determino un indice de calidad quimica del 



suelo a partir del PCA normalizado de las propiedades qulmicas del suelo 
(Velasquez et al. 2007); este Indice considero solo las variables que 
contribuyen a la fertilidad qulmica del suelo: i,Ca2 intercambiable ?, Mg2 
?, K? y NH4? y P. extralble Capacidad de almacenamiento de agua 
disponible para la planta. Mas tiempo, las mediciones de retencion de 
agua del suelo a diferentes potenciales hldricos siguieron un protocolo 
de muestreo especlfico: nucleos con estructura no perturbada (cilindros 
de 100 cm3) desde el horizonte de 0 a 10 cm (cuatro replicas) se tomaron 
de uno de los cinco puntos en cada granja, elegidos para que los 
principales tipos de uso de la tierra se muestrearan en cada area 
proporcionalmente a su ocurrencia. A partir de mediciones de laboratorio 
con un aparato de placa de presion (Pansu y Gautheyrou 2006) para 27 
puntos (108 nucleos) en cada pals, se generaron modelos de regresion 
lineal multiple (Pachepsky y Rawls 2004) para estimar, a partir de las 
variables de suelo mas simples medidas en todos los puntos (contenido de 
arcilla, limo y arena, densidades aparentes (qb), resistencia vertical 
(Rv) del horizonte superficial medido con un penetrometro de cono, 
contenido de pH, CEC y C), las capacidades de retencion de agua a 
diferentes potenciales de agua. Luego calculamos la capacidad de agua 
disponible de la planta como el volumen de agua drenado entre potenciales 
de matriz de -30 kPa y -16 MPa. Apuntamos a separar la variabilidad en el 
ES del suelo explicada por la transformacion del paisaje en campos 
cultivables de la variabilidad debido a la diversidad del suelo en el 
paisaje, que esta relacionada con la naturaleza del sustrato geologico u 
otros factores pedogeneticos. Con este fin, distinguimos las propiedades 
dinamicas o manejables del suelo versus las propiedades inherentes del 
suelo, segun lo definido por Robinson et al. (2009) y Dominati et al. 

(2010) Entre las propiedades del suelo utilizadas para calcular la ES del 
suelo Tabla 3 Metricas del servicio del ecosistema (ES) estimadas en este 
estudio y las caracteristicas del suelo o de la planta medidas para la 
estimacion de cada metrica de ES Metrica de ES Unidad de suelo y 
variables de planta C stock en biomasa de planta Mg ha-1 Sobrecarga de 
biomasa vegetal seca de arboles (BT) y arbustos (BB) C stock en el suelo 
Mg ha-1 Densidad aparente (qb) y contenido total de C a profundidades de 
0-10, 10-20 y 20-30 cm Velocidad de infiltracion de agua en el suelo mm 
h-1 Velocidad de infiltracion (Vi) Calidad quimica del suelo 0.1-1 Ca2 
intercambiable ?, Mg2 ?, K? y NH4? y contenido extraible de fosforo (P) a 
0-10 cm de profundidad. Capacidad de almacenamiento del suelo (altura) de 
agua plantable cm Contenido de arcilla, limo y arena, qb, resistencia 
vertical (Rv), contenido de C, pH, capacidad de intercambio cationico 
(CEC) a 0-10 cm de profundidad Las cursivas indican propiedades 
inherentes del suelo 318 Landscape Ecol (2014) 29: 311-328 123 metricas 
(Tabla 3), los contenidos de arcilla, limo y arena y, en menor grade, CEC 
estan relacionados con la composicion mineral de el suelo y por lo tanto 
inherentes, porque tienen poca sensibilidad a los cambios en el uso del 
suelo. En contraste, todas las otras caracteristicas, relacionadas con la 
estructura del suelo (qb, Rv, Vi) o sus composiciones organicas y 
quimicas (C, pH, A13 ?, Ca2 ?, Mg2 ?, K ?, NH4?, P) reaccionan mas 
dinamicamente a las practicas agricolas que influyen en el suelo. Al 
igual que las metricas de ES, los contenidos de arcilla, limo y arena, y 
los valores de CEC se promediaron en cada granja. Analisis estadisticos 
de los efectos del paisaje en las metricas de ES Se analizo la 
variabilidad en las metricas de ES (1) en cada punto (es decir, la escala 
de la parcela) en funcion del uso de la tierra observado en el memento de 
la medicion y (2) en la escala de la granja como una funcion de las 
clasificaciones derivadas de la composicion y estructura del mosaico de 
cobertura del suelo. Realizamos un PCA normalizado para tener en cuenta 
los gradientes que separaron los puntos y las granjas de acuerdo con las 
metricas de ES y las propiedades inherentes del suelo. Todas las metricas 
y propiedades fueron corregidas para la distribucion normal (prueba de 



normalidad de Shapiro-Wilk P [0.05). Para evaluar el grade de 
discriminacion de las categorlas de use de la tierra o cobertura de la 
tierra en las escalas de parcelas y granjas, utilizamos analisis de 
redundancia (per ejemplo, Borcard et al.l992). El Indice que caracteriza 
la separacion entre el use del suelo o los cambios en la cobertura del 
suelo (varianza entre clases) se probo con valores simulados obtenidos 
despues de 999 permutaciones de las filas de la tabla de metricas ES. 
Finalmente, realizamos un analisis de redundancia parcial para 
descomponer la variacion de las metricas de ES de acuerdo con la 
combinacion de categorlas de uso del suelo o cobertura del suelo, 
propiedades inherentes del suelo y su interaccion. Las variaciones 
resultantes se probaron de manera similar contra valores simulados 
obtenidos despues de 999 permutaciones de las filas de la tabla de 
metricas ES. En la escala de parcelas, las graficas de caja y las pruebas 
de suma de range de Kruskal-Wallis tambien se usaron para evaluar los 
efectos unicos del uso de la tierra en cada metrica de SE. La PCA 
normalizada realizada en propiedades inherentes del suelo produjo tres 
gradientes sucesivos no correlacionados (ejes de PCA) que se usaron en el 
modelado posterior. Luego realizamos el modelado lineal general gaussiano 
(GLM) en cada metrica ES utilizando el formulario: metrica ES [i] * xl? 
x2? x3? xl: zl? x2: zl? x3: zl, con xl, x2 y x3 como las primeras 
puntuaciones de la trama de tres ejes respectivamente y zl como las 
categorlas de uso de la tierra. Para la seleccion del modelo, utilizamos 
el Criterio de informacion de Akaike (AICc) y, segun lo sugerido por 
Burnham y Anderson (2002: p. 71), nos enfocamos en las diferencias en 
AICc para los primeros tres modelos. La significacion estadistica de los 
parametros del modelo se estimo mediante el analisis de la desviacion 
(validado por las pruebas F) entre un modelo nulo (sin efecto) y un 
modelo que incluia una o mas variables explicativas. El analisis de la 
desviacion tambien nos permitio calcular la proporcion de la desviacion 
explicada por las variables explicativas (es decir, un indice de la 
contribucion de la variable a la variacion metrica ES). Las estadisticas 
y los resultados graficos se calcularon con las bibliotecas ade4 (Dray y 
Dufour 2007), Mulln (Barton 2012) y veganas (Oksanen et al. 2013) 
implementadas en R freeware (R Core Team 2013) . Resultados Variacion ES 
en la escala de la parcela Los analisis estadisticos se realizaron en 260 
puntos despues de eliminar 10 puntos con valores extremes para ciertas 
variables (Tabla 2): seis puntos de Brasil con los valores mas altos de 
stock de C en biomasa vegetal (dos puntos) o suelo indice de calidad 
quimica (dos puntos) o el contenido mas bajo de arcilla (dos puntos) y 
cuatro puntos de Colombia con los valores mas altos de tasa de 
infiltracion de agua en el suelo. Los primeros dos ejes de un PCA 
realizado en las metricas ES representaron el 59% de la varianza total. 

El primer eje PCA (39% de la varianza total) contrastaba los bosques 
primarios (preservados, explotados y quemados) y los pastizales, mientras 
que el segundo eje PCA tendia a oponerse a los cultivos anuales y los 
bosques secundarios (Fig. 3a, b). La clasificacion del uso del suelo 
(Tabla 2) explico el 36% (es decir, la varianza entre clases; P \ 0.001) 
de la varianza total en ES. Los bosques primarios se asociaron con un 
mayor stock de C de biomasa vegetal y tasas de infiltracion de agua. En 
contraste, se observe una mayor capacidad de agua disponible para las 
plantas y la calidad quimica del suelo en los pastizales. El analisis de 
componentes principales realizado en la textura del suelo (arcilla, limo) 
y CEC revelo una separacion similar de los usos de la tierra a lo largo 
del primer eje (63% de la varianza total), en contraste con los bosques 
primarios y pastizales no quemados (Fig. 3c, d). La clasificacion del uso 
del suelo explico el 35% (P \ 0.001) de la variacion total de Landscape 
Ecol (2014) 29: 311-328 319 123 en las propiedades inherentes del suelo. 
Los bosques primarios se asociaron con suelos que tienen un CEC mas bajo, 
menos arcilla y poco limo. Dado que el uso de la tierra tuvo poco o 



ningun impacto en estas propiedades inherentes del suelo, PCA destaco las 
diferencias en los suelos entre los sitios, que pueden atribuirse 
principalmente a los sustratos geologicos. Los resultados de la Fig. 3 
del analisis normalizado de componentes principales (PCA) se realizaron 
en las cinco metricas de servicio del ecosistema (a, b) y las propiedades 
inherentes del suelo (c, d) a escala de parcela. a, c Correlacion de 
variables con los dos primeros ejes PCA (cada flecha apunta en la 
direccion del valor mas alto para una variable dada); El recuadro muestra 
la distribucion de la varianza entre los ejes PCA. b, d Mapa factorial 
PCA asociado que muestra parcelas agrupadas por clases de uso de la 
tierra (R2 = 0.339, P \ 0.001) y por propiedades inherentes del suelo (R2 
= 0.347, P \ 0.001) 320 Landscape Ecol (2014) 29: 311-328 123 La 
naturaleza cristalina y la edad de las rocas del Escudo Brasileno han 
favorecido la formacion de suelos arenosos en lugar de arcillosos y 
arcillosos; La fraccion arcillosa de estos suelos esta compuesta solo por 
caolinita y oxidos metalicos con bajo CEC, a diferencia de los suelos 
desarrollados en formaciones sedimentarias mas recientes del piedemonte 
andino. Como era de esperar, la variabilidad en ES resulto de una 
variedad de factores, incluidas las propiedades inherentes del suelo, el 
uso del suelo y su interaccion. 
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El analisis de redundancia parcial demostro que estos tres efectos 
explicaron el 46% de la varianza total en ES, principalmente debido al 
uso de la tierra y la interaccion entre el uso de la tierra y las 
propiedades inherentes del suelo (37% de la varianza explicable cada uno) 
(Tabla 4). El modelado lineal general realizado en cada una de las cinco 
metricas ES destaco por separado la gran contribucion de las propiedades 
inherentes del suelo a la desviacion explicada; El efecto de la 
interaccion entre ellos y el uso de la tierra fue reducido o 
insignificante (Tabla 5). Los efectos de las propiedades inherentes del 
suelo fueron mayores en la calidad qulmica e incluso mayor en la 
capacidad de agua disponible para la planta, ya que las relaciones entre 
el potencial hldrico y el contenido de agua a partir del cual se estimo 
esta ultima metrica dependen en gran medida de la textura del suelo. El 
uso de la tierra, ya sea solo o en interaccion, tuvo un impacto mucho 
menor, especialmente en la capacidad de agua disponible para las plantas 
del suelo. Un bajo efecto de las propiedades inherentes del suelo a veces 
se compenso con un mayor efecto de uso de la tierra, como para el stock 
de C en la biomasa vegetal y la tasa de infiltracion de agua (Tabla 5). 
Los diagramas de caja de las metricas de ES en funcion de las 11 
categorlas de uso de la tierra mostraron una alta variacion dentro de la 
categorla (Fig. 4). Para las existencias de C del suelo y la capacidad de 
almacenamiento de agua disponible en el suelo, las tendencias inherentes 
al uso de la tierra (prueba de Kruskal-Wallis; P \ 0.001 para todas las 
metricas de ES) a lo largo de la sucesion bosque primario-cultivos- 
pastos-plantaciones-bosque secundario fueron impulsadas por las 
propiedades inherentes del suelo (Tabla 5). En contraste, la tasa de 
infiltracion de agua, que depende del uso de la tierra, fue 
significativamente menor en las tres categorlas de pastures. Variacion de 
ES a escala de la granja El analisis del componente principal realizado 
en las metricas de ES promediadas para cada granja nuevamente revelo la 
importancia del stock de C en la biomasa de la planta, as! como la 
capacidad de agua disponible de la planta del suelo, en la separacion de 
las granjas a lo largo del primer eje (44% de varianza total; Fig. 5a). A 
diferencia del analisis a escala de parcela, la tasa de infiltracion de 



agua no vario con estas metricas. El analisis entre clases mostro que la 
clasificacion de la composicion de la cobertura del suelo en 2007 
explicaba el 35% (P \ 0.001) de la varianza total de ES. El primer eje 
PCA identified un gradiente de presidn antropogenica, que contrastaba con 
dos conjuntos de granjas, cada una de las cuales contenla tres clases que 
no diferlan mucho, a pesar de que estaban diferenciadas por el analisis 
del paisaje (Fig. 5c). Por lo tanto, el conjunto de clases Ac, Be y Cc, 
asociado con un mayor stock de C en la biomasa vegetal (Fig. 5a), inclula 
granjas brasilenas en las cuales el bosque primario todavla representaba 
45-98% del area de la granja (indicador de composicidn de 2007 entre 1 y 
0.6). En contraste, el conjunto de clases Ec, Fc y Gc mostrd una mayor 
divisidn del analisis de redundancia del suelo Tabla 4 (PDA) particidn de 
los servicios del ecosistema en la escala de la parcela de acuerdo con el 
uso del suelo, las propiedades inherentes del suelo y su interaccion 
Varianza R2% del total explicable Total explicable 2.298 *** 0.460 100 
Propiedades del suelo inherentes 0.603 *** 0.121 26 Uso del suelo 0.848 
*** 0.170 37 Uso del suelo 9 propiedades del suelo inherentes 0.847nt 
0.169 37 *** P simulado <0.001 para las pruebas de permutacion de RDA nt 
No comprobable Tabla 5 Descomposicion de desviacion obtenida de GLM 
realizada a escala de parcela entre cada metrica de servicio del 
ecosistema (ES) y categorlas de uso de la tierra y propiedades inherentes 
del suelo (los valores corresponden al porcentaje de desviacion 
explicado) Metricas de servicio del ecosistema Uso de la tierra 
Propiedades inherentes del suelo Interaccion C stock en el suelo - 22.8 
*** - Existencias de C en la biomasa vegetal 52.6 *** 12.2 *** - Tasa de 
infiltracion de agua en el suelo 22.4 *** 10.7 *** - Capacidad de 
almacenamiento de agua disponible en la planta 4.3 ** 52.8 *** 4.2 * * 
Suelo c calidad qulmica 4.2 * 41.3 *** - El GLM de cada metrica ES tiene 
la forma: metrica ES [1] * xl? x2? x3? xl: zl? x2: zl? x3: zl, con xl, x2 
y x3 como los primeros puntajes de parcela de tres ejes, respectivamente, 
del analisis normalizado de componentes principales realizado en las 
propiedades inherentes del suelo, y zl como las categorlas de uso de la 
tierra *** P <0.001; ** P <0,01; * P <0.05 Landscape Ecol (2014) 29: 311- 
328 321 123 capacidad de almacenamiento de agua disponible para la planta 
y calidad quimica (Fig. 5a). Este conjunto inclula granjas colombianas 
cuyo paisaje se vio intensamente afectado por la agricultura (indicador 
de composicion de 2007 entre 0 y 0.35), caracterizado por pasturas o 
cultivos perennes y cubiertas de tierras forestales en barbecho o 
secundario. La clase Dc, en el centre de la parcela PCA (Fig. 5c), agrupo 
granjas de areas Mac.aranduba y Aguadulce (Fig. 5b) con un paisaje de 
transicion entre las que se encuentran en cada uno de los dos conjuntos 
identificados anteriormente. La clasificacion basada en la dinamica de la 
composicion de la cubierta terrestre de 1990 a 2007 explico la mayor 
parte de la variacion de ES (45%, P \ 0.001) y distinguio tres conjuntos 
de granjas (Fig. 5d). El primer conjunto estaba compuesto por granjas de 
areas brasilenas con bosque primario restante. Las clases Acd, Bed, Ced y 
Dcd de este primer conjunto (indicador de dinamica de composicion entre 1 
y 0.5) se distinguen por la intensidad de la tala de bosques primaries 
entre 1990 y 2007. El segundo conjunto de clases. Fed y Fed (indicador de 
dinamica de composicion entre 0.2 y 0), estaba dominado por granjas 
colombianas con un paisaje agricola, donde aparecieron barbechos o 
coberturas secundarias entre 1990 y 2007. El tercer conjunto correspondio 
a algunas granjas (clase Gcd) cuyo paisaje cambio poco o nada, dominado 
por un pais homogeneo. paisaje agricola desde 1990 (pastes o cultivos 
perennes). La clasificacion basada en la estructura de cobertura del 
suelo medida en 2007 (Fig. 5e) solo explico el 23% de la varianza metrica 
total de ES (P \ 0.001) . La clasificacion de la dinamica de la estructura 
de la cobertura del suelo explico el 29% de la varianza metrica total de 
ES (P \ 0.001) y distinguio dos conjuntos de clases de granja a lo largo 
del primer eje PCA (Fig. 5f, clases que no estan en orden de indicador de 



dinamica de la estructura de la cobertura del suelo). Las clases Asd, 

Csd, Dsd y Esd se caracterizaron per una escasa fragmentacion forestal 
debido a la deforestacion reciente, o un paisaje muy heterogeneo y 
dinamico con una importancia relativa alta de la vegetacion lenosa 
(bosque primario o secundario y cultivos perennes). Las clases Bsd, Fsd y 
Gsd se caracterizaron por un paisaje agricola homogeneo y estable o una 
baja diversidad de parches en un sistema silvopastoril reciente (area de 
Balcanes). Finalmente, el analisis de redundancia parcial demostro que, 
entre las clasificaciones, la dinamica de la composicion de la cobertura 
del suelo explicaba la mayor parte de la varianza de la metrica ES. Fig. 

4 Diagramas de caja de las metricas del servicio del ecosistema en 

funcion de las 11 categorlas de uso de la tierra. Fig. 5 Resultados de un 
principle normalizado El analisis de componentes (PCA) se realizo en las 
cinco metricas de servicios del ecosistema a escala de granja. una 
correlacion de metricas ES con los dos primeros ejes PCA (cada flecha 
apunta en la direccion del valor mas alto para una metrica dada); el 
recuadro muestra la distribucion de la varianza entre los ejes PCA (es 
decir, valores propios). Mapas factoriales de PCA asociados que muestran 
granjas agrupadas de acuerdo con: b area de estudio; c composicion de la 

cobertura del suelo en 2007 (R2 = 0.354, P \ 0.001); d dinamica de la 

composicion de la cubierta terrestre desde 1990 hasta 2007 (R2 = 0.451, P 

\ 0.001); e estructura de cobertura del suelo en 2007 (R2 = 0.225, P \ 
0.001); yf dinamica de la estructura de la cubierta terrestre de 1990 a 
2007 (R2 = 0.287, P \ 0.001) c 322 Landscape Ecol (2014) 29: 311-328 123 
(45%, P \ 0.001) a escala de la granja, ambos por si mismos ( 26% de la 
varianza explicada) y en interaccion con las propiedades inherentes del 
suelo (53% de la varianza explicada) (Tabla 6). La dinamica de la 
estructura de la cobertura del suelo explico el 29% de la varianza 
metrica ES (P \ 0.01), por si misma y en interaccion con las propiedades 
inherentes del suelo Landscape Ecol (2014) 29: 311-328 323 123 (19 y 35% 
de la varianza explicada, respectivamente). En contraste, tanto la 
composicion como la estructura de la cobertura del suelo medidas en 2007 
tenian poca capacidad predictiva, esencialmente debido a la interaccion y 
a un efecto puro no significativo. El efecto puro de las propiedades 
inherentes del suelo fue significativo independientemente de la 
clasificacion de cobertura del suelo considerada y mayor con las basadas 
en la estructura (Tabla 6). Discusion Nuestro proposito era abordar un 
conjunto de ES reguladores y de apoyo de los cuales dependen los 
agricultores para administrar la tierra que han limpiado. Analizamos como 
evolucionan estos ES con cambios en el paisaje a lo largo de un gradiente 
de intensidad de uso del suelo. En primer lugar, la capacidad productiva 
del suelo, es decir, su fertilidad quimica y fisica a traves del ciclo de 
nutrientes y agua, debe mantenerse, si no mejorarse. Los suelos 
amazonicos son quimicamente pobres y fisicamente fragiles, su fertilidad 
no debe ser degradada por practicas inapropiadas para el clima y las 
condiciones del suelo. En segundo lugar, estas practicas tambien deben 
ocuparse de los servicios de regulacion del ciclo del agua para conservar 
el suelo y los recursos hidricos: la precipitacion debe infiltrarse en el 
suelo y reponer el agua disponible y el agua subterranea en lugar de 
escurrirse a la superficie del suelo y aumentar los riesgos de erosion y 
nutrientes. lixiviacion aguas arriba y de inundaciones aguas abajo. 
Finalmente, los mecanismos de compensacion, como los "creditos de 
carbono", recomendados en las negociaciones internacionales sobre el 
clima, suponen que las reservas de carbono en los compartimentos de 
agroecosistemas del suelo y la vegetacion y su variabilidad en el paisaje 
son bien conocidas (Hall 2008) . Nuestros resultados confirman la 
importante contribucion de las propiedades inherentes del suelo a la 
provision de ES del suelo y su variabilidad. La diversidad de los 
entornos biofisicos amazonicos es conocida (Hoorn et al. 2010) . Como 
resultado, los suelos difieren en escalas desde la regional (Quesada et 



al. 2012) hasta la escala local de catenas de suelo (Fritsch et al. 

2007), en las que se estudia principalmente la ecologla del paisaje. Los 
suelos difieren principalmente en su profundidad y composicion. 
Consideramos propiedades del suelo inherentes de facil acceso (textura, 
CEC) y las medimos en la capa mas superficial. Estas propiedades 
explicaron el 30% de la varianza en ES en la escala de la parcela, cuando 
se toman juntas (Tabla 4) y del 11 al 71% cuando se toman por separado 
(Tabla 5). Tal contribucion se esperaba para el ES del suelo, e incluso 
para el stock de C en la biomasa vegetal, que depende de la calidad del 
suelo (Quesada et al. 2012) . Sin embargo, la transformacion del paisaje 
tuvo un impacto significative ya que aproximadamente la mitad (46%) de la 
varianza total en ES a la escala de la parcela se explico conjuntamente 
por el use del suelo, las propiedades inherentes del suelo y su 
interaccion (Tabla 4). Todavia una proporcion considerable de la varianza 
permanecio sin explicacion, agregando la incertidumbre de la medicion que 
influye en la evaluacion de las metricas ES y los limites de las dos 
clasificaciones consideradas. De hecho, la clasificacion del uso de la 
tierra no incluyo la diversidad de las practicas agricolas o la duracion 
de los cultivos. Para las propiedades inherentes del suelo, no 
consideramos, por ejemplo, el grade de hidromorfia vinculado a la 
particion de la variacion del analisis de redundancia de la Tabla 6 (RDA) 
de las metricas del servicio del ecosistema a escala de la granja de 
acuerdo con los indicadores de cobertura del suelo y las propiedades 
inherentes del suelo Variacion R2% del total explicable Composicion de la 
cubierta terrestre en 2007 Total explicable 2.537 *** 0.507 100 
Propiedades inherentes del suelo 0.767 *** 0.153 30 Paisaje 0.409ns 0.082 
16 Paisaje 9 Propiedades inherentes del suelo 1.360nt 0.272 54 Dinamica 
de la composicion de la cubierta terrestre Total explicable 2.872 *** 
0.574 100 Propiedades del suelo inherentes 0.616 *** 0.123 21 Paisaje 
0.745 *** 0.149 26 Paisaje 9 Propiedades del suelo inherentes 1.512nt 
0.302 53 Estructura de cobertura del suelo en 2007 Total explicable 2.592 
*** 0.518 100 Propiedades del suelo inherente 1.469 *** 0.294 57 Paisaje 
0.465 ns 0.093 18 Paisaje 9 propiedades inherentes del suelo 0.659nt 
0.132 25 Dinamica de la estructura de la cubierta terrestre Total 
explicable 2.631 *** 0.526 100 Propiedades inherentes del suelo 1.195 *** 
0.239 45 Paisaje 0. 503 ** 0.101 19 Paisaje 9 propiedades inherentes del 
suelo 0.932nt 0.186 35 *** P simulado <0.001 para las pruebas de 
permutacion RDA; ** P simulado <0.01 ns No significativo (P [0.05) nt No 
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topografica. En consecuencia, el uso de la cobertura de la tierra o el 
uso de la tierra como sustitutos de ES para dibujar mapas ignora la gran 
variabilidad en las metricas de ES segun el uso de la tierra (Fig. 4). 
Aunque la variabilidad dentro del uso de la tierra fue a menudo alta, 
ocurrieron variaciones significativas entre el uso de la tierra para 
todas las metricas de ES (Fig. 4). Las metricas menos influenciadas por 
la variabilidad dentro del uso de la tierra y mas influenciadas por los 
Cambios en el uso de la tierra fueron el stock de C en la biomasa vegetal 
y la tasa de infiltracion de agua (Tabla 4). Estas metricas tambien se 
correlacionaron mas positivamente entre si y negativamente con la calidad 
quimica del suelo (Fig. 3a). Este analisis mostro claramente una sinergia 
entre el stock de C en la biomasa vegetal y la tasa de infiltracion de 
agua, y un compromise demostrado por la variacion opuesta entre ellos y 
la calidad quimica del suelo. En contraste, el stock de C en el suelo no 
se explico por el uso del suelo (Tabla 4) ni se correlaciono con el stock 
de C en la biomasa vegetal (Fig. 3a). Por lo tanto, el suelo no compensa 
significativamente la perdida de C de la biomasa vegetal (Desjardins et 
al. 2004) . Los impulsores subyacentes de la dinamica temporal de los ES 
estudiados son principalmente socioeconomicos. Los sistemas agricolas de 
los pequenos productores considerados en este estudio comprenden 
practicas forestales, agricolas y de pastoreo especificas, influenciadas 



por el acceso a los mercados, el credito y la asistencia tecnica. Estas 
practicas afectan los procesos ecologicos al modificar los habitats y las 
relaciones entre los organismos vivos (Mathieu et al. 2009). En este 
sentido, la agricultura de tala y quema cambia la calidad qulmica del 
suelo, que es baja en el bosque, ya que las cenizas de la biomasa arborea 
se incorporan al suelo. Sin embargo, es bien sabido que el reciclaje 
eficiente de nutrientes de los procesos bioticos en los ecosistemas de la 
selva tropical compensa la mala calidad qulmica de los suelos (Stark y 
Jordan 1978). Los valores mas altos se registraron asi en cultivos 
anuales, el resultado de una compensacion entre el stock de C en la 
biomasa vegetal y la calidad qulmica del suelo. A continuacion, la 
intensa actividad de los organismos del suelo en el bosque primario, en 
particular los ingenieros del ecosistema del suelo (lombrices de tierra, 
termitas y hormigas), mantiene una macroporosidad favorable a la 
infiltracion y la aireacion del suelo (Barros et al. 2001) . Bajo 
pasturas, los suelos pueden compactarse en pocos anos debido al pisoteo 
del ganado, la exposicion a la precipitacion y la disminucion de la 
biodiversidad de la fauna del suelo (Mathieu et al. 2009). La 
compactacion del suelo da como resultado una perdida de macroporosidad, 
que puede compensarse con una ganancia de mesoporosidad (es decir, poros 
medianos (0.2-15 Im) que almacenan agua disponible para las plantas), 
cuyos valores mas altos se observaron en los pastos (Fig. 3a, b ); sin 
embargo, demostramos que esta metrica de ES estaba mas influenciada por 
las propiedades inherentes del suelo que el uso del suelo (Tabla 5). Los 
Cambios en el suelo se suman a la dificultad de manejar las malezas y el 
crecimiento lenoso, explicando la disminucion en la productividad de los 
pastos y finalmente su abandons (Mitja et al. 2008) . La posicion 
intermedia de los bosques secundarios en el analisis de las metricas de 
SE, entre pasturas y bosques primaries, revela una restauracion parcial 
de las propiedades fisicas del suelo, especialmente la tasa de 
infiltracion (Fig. 3b) (Zimmermann et al. 2006) . Esta restauracion fisica 
depende de varies factores (duracion del paste, practicas de pastoreo, 
composicion y estructura del paisaje circundante) que influyen en la 
composicion de especies de las sucesiones de plantas y su capacidad para 
restaurar el suelo mediante la reactivacion de los procesos existentes en 
el bosque. La posicion intermedia de los cultivos perennes tambien indica 
que degradan el suelo menos que los pastos. Cambiar la escala no modifica 
notablemente las relaciones entre las variables, excepto el stock de C en 
la biomasa vegetal y la tasa de infiltracion de agua, que parecen menos 
correlacionadas. Promediar los cinco puntos en cada granja para realizar 
analisis a escala de granja puede no ser el enfoque mas apropiado para 
estas dos metricas, que a menudo varian de 2 a 3 ordenes de magnitud 
entre bosques y pastizales. Los paisajes del frente pionero han cambiado 
rapidamente en los ultimos 10-20 anos, y es interesante notar que, entre 
las cuatro clasificaciones a escala de granja, la que se basa en la 
dinamica de la composicion de la cubierta terrestre (1990-2007) explico 
la mayor variacion en metricas ES. La composicion de la cobertura del 
suelo realizada en 2007 en el memento de la adquisicion de dates de campo 
explicaba una menor variacion en las metricas de ES (Tabla 5). Las 
clasificaciones basadas en la composicion separaron los paisajes 
agricolas estables de los mucho mas dinamicos y diversos relacionados con 
la rapida sustitucion del ecosistema forestal por un agrosistema; 
contrastaron los dos paisajes dominantes en los dos paises. Sin embargo, 
la mayoria de las granjas que establecen sistemas agroforestales (grupo 
Gcd; Fig. 5d) se separaron de las otras granjas colombianas, 
especialmente las granjas ganaderas extensas, que implican tasas de 
infiltracion mas altas. La estructura de cobertura del suelo explicaba 
menos la varianza en las metricas de ES. Parece que las metricas que 
describen la estructura del paisaje (fragmentacion, tamano, forma, 
metricas de conectividad) influyen en la biodiversidad (por ejemplo. 



Fahrig et al. 2011) mas que la regulacion y el apoyo a Landscape Ecol 
(2014) 29: 311-328 325 123 ES estudiado aqul a traves de C secuestro en 
el suelo y la biomasa vegetal, infiltracion y almacenamiento de agua en 
el suelo y la calidad quimica del suelo, que probablemente dependen de 
las practicas agricolas aplicadas localmente (en cada parcela), en lugar 
del medio ambiente de las parcelas. Finalmente, la hipotesis de que los 
patrones de composicion de la cubierta terrestre alteran ES es cierta en 
los frentes pioneros amazonicos. Sin embargo, las metricas de ES 
evaluadas varian de diferentes maneras durante la transformacion de un 
ecosistema forestal en agrosistemas, revelando sinergias y compensaciones 
entre servicios. Los servicios de regulacion del clima, el ciclo del agua 
y la erosion del suelo son los mas afectados por la expansion de los 
pastizales en el paisaje, en detrimento de la cubierta arborea (bosques 
primaries y secundarios, cultivos perennes), por la perdida de carbono en 
la biomasa y una fuerte disminucion de la infiltracion de agua. en el 
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